Symmetry classes of alternating-sign matrices under one roof 



Greg Kuperberg| 

UC Davis 

In a previous article (gj, we derived the alternating-sign matrix (ASM) theorem from the Izergin-Korepin 
determinant |jl^, |l3| |l^] for a partition function for square ice with domain wall boundary. Here we show that 
the same argument enumerates three other symmetry classes of alternating- sign matrices: VSASMs (vertically 
symmetric ASMs), even HTSASMs (half-turn-symmetric ASMs), and even QTSASMs (quarter-turn- symmetric 
ASMs). The VSASM enumeration was conjectured by Mills; the others by Robbins [Bin. We introduce sev- 
eral new types of ASMs: UASMs (ASMs with a U-turn side), UUASMs (two U-turn sides), OSASMs (off- 
diagonally symmetric ASMs), OOSASMs (off-diagonally, off-antidiagonally symmetric), and UOSASMs (off- 
diagonally symmetric with U-turn sides). UASMs generalize VSASMs, while UUASMs generalize VHSASMs 
(vertically and horizontally symmetric ASMs) and another new class, VHPASMs (vertically and horizontally 
perverse). OSASMs, OOSASMs, and UOSASMs are related to the remaining symmetry classes of ASMs, 
namely DSASMs (diagonally symmetric), DASASMs (diagonally, anti-diagonally symmetric), and TSASMs 
(totally symmetric ASMs). We enumerate several of these new classes, and we provide several 2-enumerations 
and 3-enumerations. 

Our main technical tool is a set of multi-parameter determinant and Pfaffian formulas generalizing the Izergin- 
Korepin determinant for ASMs and the Tsuchiya determinant for UASMs [ p9| ], We evaluate specializations of 
the determinants and Pfaffians using the factor exhaustion method. 



1. INTRODUCTION 

An alternating-sign matrix (or ASM) is a matrix with en- 
tries 1, 0, and —1, such that the non-zero entries alternate in 
sign in each row and column, and such that the first and last 
non-zero entry in each row and column is 1 . Mills, Robbins, 
and Rumsey [27 \ conjectured a formula for the number of 



ASMs of order n. This formula was first proved by Zeilberger 
in 1995 [0, and later the author found a different proof Q: 

Theorem 1 (Zeilberger). There are 



A(n) 



l!4!7!---(3n-2)! 



n\(n+\)\(n + 2)\---{2n- 1)! 



nxn ASMs. 

Theorem [l] is part of a larger, unfinished structure in enu- 
merative combinatorics, much of it conjectured by Robbins 
[ fTl] , [32| | . The structure includes two types of relations between 
alternating-sign matrices and another class of combinatorial 
objects, plane partitions in boxes. (A plane partition in a box 
is an order ideal in the poset [1, . . . ,a] x [1, . . . ,b] x [1, . . . ,c]. 
They can be interpreted as a basis for the irreducible Weyl 
representation V (cX a ) of the Lie algebra sl(a + b).) One rela- 
tion is by analogy: The number of plane partitions in a given 
box is round (meaning a product of small factors, also called 
smooth), and so is the number in any given symmetry class. 
It is usually easy to conjecture an explicit product formula for 
round numbers. Likewise Robbins found that the number of 
ASMs in most (but not all!) symmetry classes also seems to 
be round. The other relation is equinumeration. Robbins also 
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found that there are the same numbers, both round and not, 
of many types of ASMs as there are other types of plane par- 
titions. To begin with there are exactly as many ASMs with 
no symmetry as there are plane partitions with full symmetry, 
and this is what Zeilberger more directly proved. 

Our proof of Theorem [l] sheds no light on plane partitions, 
but it does rely on a connection to another important struc- 
ture in quantum algebra and statistical mechanics, the Yang- 
Baxter equation. Using the Yang-Baxter equation, Izergin 
and Korepin found a determinant formula for the partition 
function of square ice with domain wall boundary conditions 
p"^ , pj| , |T^ ]. We noted that this state model is equivalent to 
certain weighted enumerations of ASMs. Although the deter- 
minant is singular at the point where all weights are equal, it 
is generically non-singular and round along a special curve of 
weights with a coordinate q such that q = 1 is the equal-weight 
point. (As defined in Section || roundness of a polynomial in 
q is stronger than smoothness.) Finally the ^-specialization 
of the Izergin-Korepin determinant independently generalizes 
to a round determinant in two parameters p and q. The two- 
parameter determinant can be evaluated by factor exhaustion 
in p. 

In this article we generalize this argument to some of the 
previously known classes of ASMs and also some new ones. 
The Izergin-Korepin determinant may have seemed acciden- 
tal, but we find a similar formula for each symmetry class of 
ASMs after modifying the boundary conditions at the sym- 
metry lines (Theorem |l0|). (One case was found previously 
by Tsuchiya |]3^].) The round ^-specialization and its two- 
parameter generalization may also have seemed accidental, 
but we find many such specializations complemented by four 
two-parameter determinants and two three-parameter Pfaffi- 
ans (Sections || and 0). Besides ordinary enumeration, the q- 
specializations also include several jc-enumerations (in which 
each orbit of — 1 entries in the ASM has weight x) with x = 2 
and x = 3. (Along the way we will correct an error in the 3- 
enumeration of ASMs in Reference O found by Robin Chap- 
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man.) Besides the ^-specializations we establish that many of 
the x-enumerations divide each other or otherwise share large 
factors (Section [7]). 

We speculate that our constructions are part of a yet larger 
and possibly less exotic structure in quantum algebra. In par- 
ticular the solution to the Yang-Baxter equation that we use 
corresponds to the 2-dimensional representation of the Lie al- 
gebra sl(2). We have not investigated what happens when 
sl(2) is replaced by another Lie algebra or the 2-dimensional 
representation by another representation. 
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2. STATEMENT OF RESULTS 

In general the x- enumeration of a class of ASMs is de- 
fined as the total x-weight. The x-weight of an ASM is x" 
if it has n symmetry orbits of negative entries. Figure [l] 
shows an example. (In the figures, we use + and — for 1 
and —1.) Note that the 2-enumerations of any class of ASMs 
is a special case (called the free fermion point in statistical 
mechanics) that can often be established byother methods 
% |, H 0, [l| H H ||, HQ]. The reader 

can also consider the elementary case of O-enumeration since 
it is also often round. 
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Figure 1 : A VS ASM with x-weig htx 2 . 



We summarize the conventional symmetry classes of ASMs 
and what is known and conjectured about their basic enumer- 
ations. 



• ASMs - alternating-sign matrices. The 1-, 2-, and 3- 

enumerations are all previously known. 

• VSASMs - vertically symmetric ASMs. The 2-enumeration 

is previously known. Mills conjectured the 1- 

enumeration [fffi]. We establish the 1-, 2-, and 3- 
enumerations. 

• HTSASMs - half-turn symmetric ASMs. Robbins con- 

jectured the 1 -enumeration and established the 2- 

enumeration [Mh. We establish the 1-, 2-, and 3- 
enumerations, but only for even order. 

• QTASMs - quarter-turn symmetric ASMs. Robbins con- 

jectured the enumeration. We establish the 1- and 2- 
enumerations only for even order. We also prove a for- 
mula for a factor of the 3-enumeration; the other factor 
does not appear to be round. 

• VHSASMs - vertically and horizontally symmetric ASMs. 

Robbins conjectured the 1 -enumeration and the 2- 
enumeration is previously known [O]. We only estab- 
lish a determinant formula. 

• DSASMs - diagonally symmetric ASMs. Their number 

does not appear to be round and no determinant or Pfaf- 
fian formula is known. 

• DASASMs - diagonally and antidiagonally symmetric 

ASMs. Robbins conjectured the enumeration for odd 
order, but their number does not appear to be round 
for even order. No determinant or Pfaffian formula is 
known. 

• TSASMs - totally symmetric ASMs. Their number does 

not appear to be round and no determinant or Pfaffian 
formula is known. 

Our results for these classes are as follows: 
Theorem 2. The number ofnxn ASMs is given by 

A(»)=(-3)© n ^=^. 

i.f J-i+n 

The number of 2n + 1 x 2n + 1 vertically symmetric ASMs 
(VSASMs) is given by 

A v (2n + l) = (-3f II 3U ~Jl + ' 

i,}<2n+lJ- l + 2n+1 

The number of 2n x 2n half-turn symmetric ASMs (HTSASMs) 
is given by 

a ht( 2 ») = r r>(*'TT 3 k'~'') +2 

A{n) 1 V [j j-i + n ■ 

The number of An x An quarter-turn symmetric ASMs ( QT- 
SASMs) is given by 



A QT (An)=A HT (2n)A(n) 2 . 
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In the statement of Theorem || and throughout this article, 
subscripts and products range from 1 to n unless otherwise 
specified. 

Theorem 3. The 2- and 3 -enumerations of ASMs and 
VSASMs are given by 



A{n;2) = 2© 
A(n;3) 



3" -" j-j 3(7-i) + l 



2 „2_„ ii 
2\j-i 



A v (2n + 1;2) = 2" ~ n 



3 2»- 



A v (2n + 1;3) = 



3(7-i) + l 
2 2n2+n i,M + i 3(7-0 

2fi,2|j 



n 



77ze 2 -enumerations of even-sized HTSASMs and QTSASMs 
are given by 



A HT (2«;2,1)=2" 2 1~[ 



2(7-0 + 1 



ij 2C/-0 

A OT (4»;2) = (-l)G)2 2 " 2 -"n 4 °'~ + 1 - 



integer n. It satisfies the alternating-sign condition and 
it has the same symmetries as a VHSASM, except that 
the central entry (*) has the opposite sign when read 
horizontally as when read vertically. The simplest VH- 
PASM is given in Figure [| 

• UASMs -2nx2n ASMs with U-turn boundary on the right. 
Figure [3] shows an example of a UASM. As the exam- 
ple indicates, a UASM is vertically just like an ASM. 
Horizontally the signs alternate if we read the 2k — 1st 
row from left to right, and then continue to alternate 
if we read the 2&th row from right to left. UASMs 
were first considered by Tsuchiya [p9|]. They generalize 
VSASMs. 



UUASMs 2 -2nx2n ASMs with U-turn boundary on the top 
and right. UUASMs generalize both VHSASMs and 
VHPASMs. 



• OSASMs - off-diagonally symmetric ASMs. I.e., DSASMs 

with a null diagonal. 

• OOSASMs - off-diagonally, off-antidiagonally symmetric 

ASMs. I.e., DAS ASMs with null diagonals. 



/0 + 0\ 
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Figure 2: The simplest VHPASM. 
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Figure 3: A UASM. 



We will also consider the following new types of ASMs: 

• VHPASMs 1 - vertically and horizontally perverse ASMs. 
A VHPASM has dimensions 4n + 1 x 4n + 3 for some 



• UOSASMs - off-diagonally symmetric UASMs. They in- 
clude TSASMs with null diagonals. 



Our remaining results involve weighted enumerations that 
are more general than just x-enumeration. We define the y- 
weight of a 2n x 2n HTSASM to be y k if the HTSASM has 
k non-zero entries in the upper left quadrant. This yields the 
(x, v)-enumeration of HTSASMs, which is round when y is 
— 1 and x is 1 or 3.) We define the x-weight of a UASM or 
a UUASM be the number of —Is, as before. We define the 
y-weight of a UASM to be y k if k of the U-turns are oriented 
upward in the corresponding square ice state. We define the 
y-weight of a UUASM the same way using the U-turns on the 
right, and define the z-weig ht of a UUASM to be z k if k of the 
U-turns on the top are oriented to the right. Thus we can con- 
sider the (x,y)-enumeration of UASMs and the (x,y,z) -weight 
of UUASMs. As with UASMs, the y-weight of a UOSASM 
is y k if k of the U-turns on the top are oriented to the right. 
By contrast, the y-weight of an OOSASM is y k if there are 2k 
more 1 s than — 1 s in the upper left quadrant. Finally we index 
the generating function of a given type of ASM by the length 
of one of its rows, counting the length twice if the row takes a 
U-turn, and we include x-, y-, and z-weight where applicable. 
For example A uo (8n;x, y) is the weighted number of An x An 
UOSASMs. 



1 Also known as j8-ASMs, since their boundary conditions are incompatible 
with VHS. 



Also known as Unix-to-Unix ASMs. 
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Theorem 4. There exist polynomials satisfying the equations 



A(2n;x) 


= 2A v (2n + 1 ; x) A v (2n ; x) 


A{2n + \;x) 


= A v (2n + l;x)A v (2n + 2;x) 


A\j(2n;x,y) 


= (y + l)"A v (2n + l;jc) 


A HT (2«;x, ±1) 


= A{n;x)A^{2n;x,±l) 


A vv (4n;x,y,z) 


= A v (2n + l;x)AW(4n;x,y,z) 


A 0Q (4n;x,y) 


= A (2n;x)AW(4n;x,y) 




= AW(4«;x)Ag(4n;x) 


A VQ (Sn;x,y) 


= AW(8«;x)A[ J 2 )(8«;x,y) 


A<&(4n;x,l) 


= A^(4n;*,l,l)A«(4»;x) 


A^(4n + 2;x,l) 


= 2A[ 2 )(4«;x,l,l)A[ 2 )(4»+4;x) 


A®(4n;x,-1) 


= (-x)"AW (4n;xf 


Ag)(8«;x,-1) 


= (-x)"A[ J 1 )(8«;x,l)A[ J 1 )(8«;x,l) 



edges enter and leave every tetravalent vertex. In particular if 
G is locally a square grid, then the set of ice states is called 
square ice [|25||. More generally G may also have some uni- 
valent vertices, which are called boundary, and restrictions on 
the orientations of the boundary edges are called boundary 
conditions. 



y\ yi y-i y 4 



Figure 4: Square ice with domain wall boundary. 



Many of the factorizations in Theorem Q were conjec- 
tured experimentally by David Robbins [pip; the formula for 
A v (2n;x,y) was conjectured by Cohn and Propp [||. 

Theorem 5. The generating functions in Theorem ^have the 
following special values. 



A (2«)=A V (2« + 1) 



A<$J4n; 1,1,1) 



3(7-0+2 



A^(4n;2,l,l) 



A( 2 l (4n + 2;1) 



VHP V 



= (-3)- 2 2 - n , 

i.j<2n+lJ-'+ 2n+1 

2|; 

2 „(n+2) FT 2 0'~'') + 1 

'" y4An 20' -0 

2fi,2|; 

= A V (2« + 1) 



A«(4«)=A(») 2 



AW(4«;2)-(-l)(")2"("- 1 )n 



Ag(4n;3) 



4Q-0 + 1 
j-i + n 



3&A(n) 



A«(8n)=A v (2n+l) 



AP)(8»)=A uu (4n). 

Other identities, for example that 

Ag(4n) =A HT (2n), 

are implied by combining Theorems ||, 0, Q and |J although 
such combinations do not always reflect the logic of the 
proofs. 
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Figure 5: Replacing square ice with alternating-sign entries. 



For example, a finite square region of square ice can have 
domain wall boundary, defined as in at the sides and out at the 
top and bottom, as in Figure |j. These boundary conditions 
were first considered by Korepin [ [l2| [Ts] , [l^ ]. A square ice 
state on this region yields a matrix if we replace each vertex 
by a number according Figure |[ It is easy to check that this 
transformation is a bijection between square ice with domain 
wall boundary and alternating-sign matrices QKl [23|]. 



Figure 6: Square ice with VS boundary. 



3. SQUARE ICE 

If G is a tetravalent graph, an ice state (also called a six- 
vertex state) of G is an orientation of the edges such that two 



There are also easy Injections from ice states of the graphs 
in Figures g, §, and |f to the sets of VS ASMs, VHS ASMs, 
even HTSASMs, and even QTSASMs. (The labels in these 
figures will be used later.) The dashed line in the QTSASM 
graph means that the orientation of an edge reverses as it 
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±x, 



±x 




Figure 7: Square ice with HTS boundary. 



Figure 8: Square ice with VHS boundary. 
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y\ yi 

Figure 10: Square ice with U boundary. 
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y\ y x yi y 2 



Figure 11: Square ice with UU boundary. 



crosses the line. The HTSASM and QTSASM graphs are ob- 
tained by quotienting the unrestricted ASM graph by the sym- 
metry. The median of a 2n + 1 x In + 1 VSASM is always 
the same, so we can delete it and consider the alternating-sign 
patterns on the left half. The deleted median then produces 
the alternating boundary in Figure |[ Likewise we can quarter 
a VHSASM by deleting both medians, which produces two 
alternating sides. 

Finally the square ice grids corresponding to UASMs, 
UUASMs, OSASMs, OOSASMs, and UOSASMs are shown 




1 3 



in Figures [k], [IT], [T2|, [l3|, and [l4]. The last three grids have 
right-angled divalent vertices; we require the orientations of a 
square ice state to either be both in or both out at these ver- 
tices. In contrast at the U-turn vertices one edge must point in 
and one must point out. 

4. LOCAL CONCERNS 

Throughout the article we assume the following abbrevia- 
tions: 

x =x~ l 
o(x) =x — X 
cc(x) = a (ax) a (ax). 



H 

x 2 



Figure 9: Square ice with QTS boundary. 



Figure 12: Square ice with OS boundary. 
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Figure 13: Square ice with OOS boundary. 
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Figure 15: Weights for vertices and U-turns. 



tations of internal edges. For example, the graph 



Figure 14: Square ice with UOS boundary. 



(As we discuss below, a is a global parameter that need not 
appear as an explicit argument of a.) 

We will consider a class of multiplicative weights for sym- 
metric ASMs. By a multiplicative weight we mean that the 
weight of some object is the product of the weights of its 
parts. In statistical mechanics, multiplicative weights are 
called Boltzmann weights, and the total weight of all objects 
is called a partition function. Figure |l5| shows the weights 
that we will use for the six possible states of a vertex. The 
figure also shows the weights for U-turns and corners that 
are labelled with a dot; bare edges and curves have the triv- 
ial weight 1 . The vertex weights are called an R-matrix and 
the U-turn and corner weights are called K-matrices. Vertex 
and U-turn weights depend on a parameter x (the spectral pa- 
rameter) which may be different for different vertices or U- 
turns, so we will label sites by the value of x. The weights 
also depend on three parameters a, b, and c which will be the 
same for all elements of any single square ice grid, so these 
parameters do not appear as labels. 

We will use a graph with labelled vertices as a notation for 
its corresponding partition function. If the graph has unori- 
ented boundary edges, then the partition function is also inter- 
preted as a function of the orientations of the edges. On the 
other hand, our definitions imply that we sum over the orien- 



denotes the following function on the set of four orientations 
of the boundary: 



xd >o yo 







<j(a) + o(ax) a(a 2 ) + a(ax) 







In this notation a vertex is not quite invariant under rota- 
tion by 90 degrees, so the meaning of a label depends on the 
quadrant in which it appears. The following relation holds: 



As a further abbreviation, if we label two lines of a graph 
that cross at an unlabelled vertex, the spectral parameter is set 
to their ratio: 



xy 



The labelled graphs in Figures § §, 0, g, |o[ [Tl], |T^, g and 
[l4] then represent the partition functions 

Z(n;x,y) Z^ T (n;xJ) Z v (n;x,y) Z vv (n;x,y) 
ZQj(n;x) Z (n;x) Z 00 (n;x) Z u0 (n;x) 
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Here the vectors x and y have length n when both are present, 
and otherwise x has length 2n. In the HT and OO cases there 
is a single extra parameter taken from the set {+,—}; if it is 
— then the spectral parameters in the upper half of the grid are 
negated. (Note that the index n is not defined in the same way 
as for the enumerators such as A m (2n).) 

The key property of the 7?-matrix is that it satisfies the Yang- 
Baxter equation: 

Lemma 6 (Yang-Baxter equation). Ifxyz = a, then 




As usual the Yang-Baxter equation appears to be a massive 
coincidence. In our previous review of the Yang-Baxter equa- 
tion the 7?-matrix was normalized to have a particular 
symmetry: It was the matrix of an invariant tensor over the 2- 
dimensional representation of the quantum group U q (sl(2)), 
with q related to our present parameter a. This symmetry re- 
duced the coincidence in the equation to a single numerical 
equality. The spectral parameters were chosen to satisfy the 
equality. Here we normalize the /^-matrix to reveal combina- 
torial symmetry rather than symmetry from quantum algebra. 

Proof. Taken literally, the equation consists of 64 numerical 
equalities, because there are 64 ways to orient the six bound- 
ary edges on each side. However, both sides are zero unless 
three edges point in and three point out. This leaves 20 non- 
zero equations. The equation also has three kinds of symme- 
try: The right side is the left side rotated by 180 degrees, all 
arrows may be reversed, and both sides may be rotated by 120 
degrees if the variables x, y, and z are cyclically permuted. 
By the three symmetries, 8 of the non-zero equations are tau- 
tological, and the other 12 are all equivalent. One of the 12 
non-trivial equations is 




In algebraic form, the equation is 

O (ay) O '{a 2 ) a (ax) — o(az)<j(a 2 ) 2 + a (ax) (7 (ay) a (a 2 ). 

Cancelling a factor of <r(a 2 ), expanding, and cancelling terms 
yields 

2-- -2 3 - - -3- 2 -2-- 

a xy + a xy = a z — za — za + a z + a xy + a xy, 

which is implied by the condition xyz = a. □ 

We will need the reflection equation [|3|, ^, [35p, an analogue 
of the Yang-Baxter equation that relates a Zf-matrix to the R- 
matrix. 



Lemma 7 (Reflection equation). If st = ay and st = ax, then 




Proof. The argument is similar to that for Lemma [| Both 
sides are zero unless two boundary edges point in and two 
point out. There is a symmetry exchanging the two sides given 
by reflecting through a horizontal line and simultaneously re- 
versing all arrows. (Note that the weights of a U-turn are not 
invariant under reflection alone.) Under this symmetry 4 of 
the 6 non-zero equations are tautological, and the other 2 are 
equivalent. One of these is: 




Algebraically, the equation reads: 

a(bt)a(a 2 )(o(ay)o(bs) + o(ax)a(bs)) = 

a (bt) a (a 2 ) (a (ax) o (bs) + a(ay)o(bs)). 

All terms of the equation match or cancel when st = ay and 
st = ax. □ 

The corner /f-matrix also satisfies the reflection equation 
§• 

Lemma 8. For any x and y, 




Proof. Diagonal reflection exchanges the two sides. Both 
sides are zero if an odd number of boundary edges point in- 
ward. If two boundary edges point in and the other two point 
out, then arrow reversal is also a symmetry, because one cor- 
ner must have inward arrows and the other outward arrows. 
These facts together imply that all cases of the equation are 
null or tautological. □ 

Finally we will need an equation that, loosely speaking, in- 
verts a U-turn: 

Lemma 9 (Fish equation). For any a and x, 

ax = a(a 2 x 2 ) ~y ax 

The proof is elementary. 
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5. DETERMINANTS 

In this section we will establish determinant and Pfaffian 
formulas for the partition functions defined in Sections || and 

Recall that the Pfaffian of an antisymmetric 2n x 2n matrix 
A is defined as 

7CEX i 

where X C S 2n has one representative in each coset of the 
wreath product S 2 lS„. (Thus X admits a bijection with the 
set of perfect matchings of {1, . . . ,2n}.) Recall also that 

detA = (PfA) 2 . 



Theorem 10. Let 



±- 



M(n;x,y) t j 

M„ T (n;x,y)- ■ ■ _ 
'" >hJ o{ax iyj ) o(axff } ) 

. _ 1 1 

M u («;x,y). . = — — — - — — — 

a \ x ffj) a \ x iyj) 

Mvufcxjjjj = _ _ _ 

o(b yj )a(c Xi ) a(foy,)tt(c^) 

i 1 L— 



M (n;x)u 



y ' a{x iXj ) 
a(x iXj ) 



c 2 b 2 



M 00 («;x),. ■ =a{x i x 1 ) 

J ' \o(aXjXj) a{aXfXj) t 

M\} ) Jn;x) l , =a{x i x.)a{x i x j ) ( -|- 

/ a{cXj)a{cXf) a(cxi)a{cx,) 
Mg(n;x) iJ =a(x i x J )a(xpc j )l ' ' 



a{ax i Xj) a{ax i Xj) 
c (cx i ) a (cxj ) a (cx { ) (7 (cxj ) 1 
a(aXjXj) <j(aXjXj) 



Then 



C7(a 2 )"n,/«(^,) 



n i <jff(¥>(}'ij'j) 
Z ht("^) - n ' J -\ 2 (detM)(detM^ T ) 



u{ n;x,y) - ^ . )n . ; a^a^ffj) 

■ (detMu) 



z uu (";^3 ? ) 



a(a 2 )"U i (y(a 2 x 2 )(y(a 2 y 2 )U l , j a(x i y j ) 2 a(x i y j ) 2 

Ui<j o^pcjYaiy^jY Ui<j a^x^aiy^j) 2 
•(detM u )(detM uu ) 

g( a 2 )»g( a ) 3 "n, <7 < 2n «(^;) 2 
n, <> <2„^(^ i ) 2 

•(PfMg)(PfMg) 



^(a 2 )" n,v,<2n 

z ° (n;f) = n aL\ (PfMo) 

Z Oo("^) = rr S,A2 (PfM Q )(PfM 00 ) 



Z uo («;x) 



n, <i <2„^(^ ; ) 2 

fc 2 "g( a 2 )"g( a ) 2 "nK 2n C7(a 2 ^ 2 ) 

n (<; < 2 „ ^(^) 2 n i < i < 2 „ °( x i x j) 2 

• n a(i,x 7 ) 2 a(x,x / ) 2 (PfM( J 1 ))(PfM( 2 ) 

i< j<2n 



We call the first four partition functions the determinant 
partition functions and the other four the Pfaffian partition 
functions. 

Remark. The partition function Z (J (n;5 ? ,y), the Tsuchiya de- 
terminant, is nearly invariant if x is exchanged with y. Sim- 
ilarly Z Q (n;x) is nearly invariant if each x i is replaced with 
x t . We have no direct explanation for these symmetries. Note 
that the first symmetry is less apparent in Tsuchiya's matrix 
M [[?9| eq. (42)], which has an asymmetric factor 

sinh(C_+A / ) sinh(C_-A ; ) 



sinh(Aj 



CO; 



sinh(Aj 



co- 



in this expression (Op A-, and £_ are obtained fromx ( , y., and 
b by reparameterization. If we factor this expression, 

sinh(2A ■) sinh(£_ - <a,) 

F- ■ = 

lJ sinh(A- + (» 1 -)sinh(A. — 0) ; ) 

we can then pull the asymmetric factors out of the determinant 
since they each depend on only one of the two indices i and j. 
This also explains why the /f-matrix parameter or b need 
not appear in the matrix M v . 



The proof of Theorem 10 uses recurrence relations that de- 
termine both sides. The relations are expressed in Lemmas |llj 
|l2| , pj| , and [l4]. Indeed, the first three of these lemmas are ob- 
vious for the right-hand sides of Theorem HQ; only Lemma [l4| 
needs to be argued for both sides. 

Lemma 11 (Baxter, Sklyanin). Each of the partition func- 
tions in Theorem 10 is symmetric in the coordinates ofx. Each 
determinant partition function is symmetric in the coordinates 
ofy. The partition functions Z v (n;x,y) and Z uu (n;.x,;y) gain 
afactorofo(a 2 x 2 )/o{a 2 x 2 ) ifx { is replaced by x i for a single 
i. Similarly Z^(n;x,y) gains (j(a 2 yj) / <r(a 2 yf) under y i t— > y i 
andZ vo {n;x) gains o(a 2 x 2 ) / o{a 2 x 2 ) underx i i— » jc f . 
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Proof. Invariance of Z(n;x,y) is an illustrative case. We ex- 
change Xj with x i+l for any i < n — 1 by crossing the corre- 
sponding lines at the left side. If the spectral parameter of the 
crossing is z = ax j x i+l , we can move it to the right side using 
the Yang-Baxter equation (Lemma pb and then remove it: 



Oiaz 



v i+l 



The argument for symmetry in x is exactly the same for all of 
the square ice grids without U-turns. If the grid has diagonal 
boundary with corner vertices, we can bounce the crossing off 
of it using Lemma [| 

If the grid has U-turn boundary on the right, we exchange 
x t with x j+l by crossing the x t line over the two lines above 
it. We let the spectral parameters of these two crossings be 
Z = dx~jX i+l and w = ax j x j+l . We move both crossings to the 
right using the Yang-Baxter equation, then we bounce them 
off of the U-turns using the reflection equation (Lemma [7]): 




Also if the grid has a U-turn on the right, we establish covari- 
ance under X: 1— > x, by switching the lines with these two labels 
and eating the crossing using the fish equation (Lemma fk. 



The same arguments establish symmetry in the coordinates 
of y. All of the arguments used in combination establish the 
claimed properties of Z uo (n;x). □ 

Lemma 12. The partition function Z HT (n;x,y) gains a factor 
o/(±l)" ifx i andy i are replaced by x, and y j for all i simulta- 
neously. Similarly Z QO ± {n;x) is invariant under x i 1— > x i and 
b 1— * c 1— » b. 

Proof. In both cases, the symmetry is effected by reflecting 
the square ice grid or the alternating-sign matrices through a 
horizontal line. □ 

For a vectors = (x i , . . . ,x n ), let x 1 = (jc 2 , . . . ,x n ). 

Lemma 13. lfx x = ay v then 



Z(n;x,y) 



Z(n-l;Z,f 
^ ( 



a{a 2 )Y\(y{ax i y x )a{ax l y i ) 



2<i 



Z± T (n;x,y) 



Z v (n;x,y) 
Z v (n- l;x',f) 

Z m (n;x,y) 
Z^in-l-IJ) 



2<i 



a(ar)a{arxi)a{by{) 

■ n a ( fli ^ 1 ) ° ( a * i^i ) ° ( ax iy 1 ) a 1^1 ) 



2<i 



a(a ) a{a x^)a{a y{)a{by^)a{cX]) 



2<i 



■Waiax^faiax-Ji) 

2<i 



<j(a) 2 <j(a 2 ) 2 Y\ <y (ax jX^ 2 0(0x^2 

3<i<2n 



- o(a 2 ) Y\ <5{ax l x i )o(ax2X l ) 

3<i<2« 

: c 2 <7(a 2 ) 2 Y\ (y(ax l x j ) 2 G(ax 2 x i ) 2 

3</<2h 

: b 2 a{a) 2 a{a 2 ) 2 a{a 2 x\)a{a 2 x\) 

• o{cX\)o(cx2) Y\ <y{ax l x i ) 2 o{ax2X i ) 

3<<<2n 

• Y\ (y{ax l x i ) 2 o{ax2X i ) 2 . 

3<K2n 



Proof. This lemma is clearer in the alternating-sign matrix 
model than it is in the square ice model. The partition func- 
tion Z(n;x,y) is a sum over n xn alternating-sign matrices in 
which each entry of the matrix has a multiplicative weight. 
When yj = ax l7 the weight of a in the southwest corner is 
0. Consequently this corner is forced to be 1 and the left col- 
umn and bottom row are forced to be 0, as in Figure The 



Ifx 2 = ax x , then 

Z/^yY ^ TT f JC ^ 

Z QT (n- l;x") 
Ifx 2 = ax x , then 

Z (n;x) 
Z (n-l;x") 

z oo(»;*) 

Z 00 (n-1;F) 

Z UQ (»;x) 
Z uo («-l;x") 
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Figure 16: ASM entries forced by y l = ax^. 



Function 


Width 


Specializations 


Z(n;x,y) 


n-\ 


y\ 


= ax i 


Z± T (n;x,y) 


2n-\ 


y\ 


= a Xf 


Z v (n;x,y) 


2n-l 


y\ 


= ax^ 




An 


y\ 


= a ±l x ±1 ,a 


Z^-^y Yl , JC ^ 


An -3 


x i 


= a Xf 


Z Q (n;x) 


In -2 


x \ 


= aXj 


Z 0Q (n;x) 


An -3 


H 


= a Xf 


Z VQ (n\x) 


8n-4 


x \ 


= a ±1 * ±1 ,a 



sum reduces to one over (n — 1) x (n — 1) ASMs. The only 
discrepancy between Z(n;x*,y)| Ji=aX] and Z(n~ l;x* ,y*) is the 
weights of the forced entries, which the lemma lists as factors. 

The argument in the other determinant cases is identical. 
The argument in the Pfaffian cases is only slightly different: 
All QTSASMs have zeroes in the corners, and the specializa- 
tion x 2 = ax l instead forces a 1 next to each corner and zeroes 
the first two rows and columns from each edge. Likewise the 
specialization x 2 = ax x forces a 1 next to each corner of an 
OSASM or an OOSASM and a 1 in the third row entry bot- 
tom of a UOSASM, and several rows and columns of zeroes 
in each of these cases. □ 

Define the width of a Laurent polynomial to be the differ- 
ence in degree between the leading and trailing terms. (For 
example, q 3 — q~ 2 has width 5.) 

Lemma 14. Both sides of each equation of Theorem [/^ are 
Laurent polynomials in each coordinate ofx ( and y in the de- 
terminant cases) and their widths in x x (y x in the determinant 
cases) are as given in Table [7J 

To conclude the proof of Theorem |l^, we claim that Lem- 
mas 11, 12|, and [j~3| inductively determine both sides by La- 
grange interpolation. (To begin the induction each partition 
function is set to 1 when n = 0.) If a Laurent polynomial of 
width w has prespecified leading and trailing exponents, it is 
determined by w + 1 of its values. Each of our partition func- 
tions is a centered Laurent polynomial in Xj (in the Pfaffian 
cases) or y l (in the determinant cases). Moreover each is ei- 
ther an even function or an odd function. Thus we only need 
w + 1 specializations, where w is the width in x\ (or y\). 

These widths are summarized in Table [jj. To compute them, 
observe that each entry in the bottom row of an ASM con- 
tributes 1 to the width. In the UASM and UUASM cases, it 
is the bottom two rows, and the U-turn itself contributes 1 
to the width as well. In the QTSASM case, the corner en- 
tries always have weight a(a) and do not contribute to the 
width. Lemmas [11] and [l3| together provide many specializa- 
tions which are listed in Table [j]. Note that Lemma [Tl| implies 
that a{a 2 x\) divides Z uu («;x,y) Z u0 (n;x, vy), which provides 
an extra specialization in these two cases. In conclusion, it is 
easy to check that there are enough specializations to match 
the widths. 

Remark. The formulas in Theorem [Io| are even more special 
than Lemmas [IT] through |l4| suggest. Among the evidence for 
this, the recurrence relations still hold with only slight modifi- 
cations if all spectral parameters in the QT, UU, and UO grids 



Table 1: Widths and specializations of partition functions.. 



are multiplied by an extra parameter z. Similarly the spectral 
parameters in the top halves of the HT and OO grids may be 
multiplied by an arbitrary z instead of by ± 1 . However, we 
were not able to generalize Theorem |l0| to include this param- 
eter. 

Lemma [l3] reveals another subtlety, namely that 

Z+ T (n;x,y) _ ( Z{n;x,y) 
Z+ T («-l;x',/) ~ \Z{n-\;x',f) 

at every specialization y i — a ±l Xj. Since this coincidence 
holds for enough specializations to determine Z^ T (n;x,y) en- 
tirely, one might suppose that 

Zft T (n;x,y) — Z(n;x,y) 2 . 

But then Z^ T (n;x,y) would be an even function of y lt while 
in reality it is an odd function. The other symmetry classes 
involving half-turn rotation have similar behavior. 

6. FACTOR EXHAUSTION 

In this section we derive several round determinants and 
Pfaffians depending on two and three parameters. We will 
later identify special cases of the determinants and Pfaffians 
with those appearing in Theorem [l(], and they will specialize 
further to establish the enumerations in Theorems |[ |, and ||. 

Since the formulas in this section may seem complicated, 
we recommend verifying that they are round without worrying 
about their exact form in the first reading. For this purpose we 
give a more precise definition of roundness that also applies 
to polynomials. A term R„ in a sequence of rational polyno- 
mials depending on one or more variables is round if it is a 
ratio of products of constants, monomials, and differences of 
two monic terms. All exponents and constant factors should 
grow poly normally in n or be independent of n. For exam- 
ple, n\3"(q + p") is round and Gaussian binomial coefficients 
are round. Note that a round polynomial in a single variable 
must be a product of cyclotomic polynomials, which is part of 
the motivation for the term "round". Roundness is preserved 
when a variable is set to 1 or to a product of other variables. 
In a later reading one can verify the explicit formulas. This is 
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a tedious but elementary computation, because all round ex- 
pressions involved have an explicit and regular form. As a 
warmup the reader can verify that the expressions for A(n) in 
Theorems [j] and [| coincide. 

We begin with the classic Cauchy double alternant and a 
Pfaffian generalization found independe ntly by Stembridge 
and by Laksov, Lascoux, and Thorup [|l6| |24 37 38 1 . 



Theorem 16. Let 



Theorem 15 (Cauchy, S. L. L. T.). Let 



1 



x i +y j 



For 1 < z, j < In, let 



C 4 (*),',; : 



Xt+Xj 



l-xpcj 



Then 



detC, 



detC, = 



PfC, 



PfC A 



Ui<j ( 1 - xfCj) ( 1 - ypj) (xj - x t ) (jj - y .) 

X\i,Mi+yj)(X+^j) 
■W+^+yi) 

n 

i<j<2n X i+ X j 



n r 

;</<2« 1 



Proof. Our proof is by the factor exhaustion method [pop. The 
determinant detC : is divisible by x= —x i because when x i = x,-, 
two rows of Cj are proportional. Likewise it is also divisible 
by y j — y { . At the same time, the polynomial 

i-j 

has degree n 2 — n, so it has no room for other non-constant 
factors. This determines detCj up to a constant, which can be 
found inductively by setting x l = —y v 

The determinant detC 2 is argued the same way. The Pfaf- 
fians PfC 3 and PfC 4 are also argued the same way; here the 
constant factor can be found by setting x 1 =x 2 - □ 

Next we evaluate four determinants in the variables p and a. 
We use two more functions similar to a and a from Sections: 



Y(q)=q l/2 -q- l/2 



■q 



1/2. 



-q 



-1/2 



rt \ 

7 { q » + j +i ) r(q n+ n 

3\P'4)iJ y^pn+j+i} y(pn+j-i) 

wqhj - T{pj+i) T{p 



J-l) 



Then 



detr, 



n ¥j r(p J •ni^qp 1 ■> 



u U 7(p n+j -') 

2 n U r(p j - l ) 2 U y{qp ] -') 

detr 2 = (-1)© 

H<;<2« 7(/> M ) LX-j^n+l 7(^ ; "0 
2L/ 

2"n i<; <„ rip^-^n,,^ 7(qp J -') 



det7; 



detr, 



I\ L j<p J - , W + ') 



Proof. Factor exhaustion. We first view each determinant as a 
fractional Laurent polynomial in q. By choosing special val- 
ues of q, we will find enough factors in each determinant to 
account for their entire width, thus determining them up to 
a rational factor R(p). (Each determinant is a centered Lau- 
rent polynomial in q with fractional exponents. The notion of 
width make sense for these.) We will derive this factor by a 
separate method. 

For example, if < k < n, then detTj is divisible by 
y(qp~ k )"~ k because 



! -<e<t 



Evidently T 1 (p,p k ) is a sum of k rank 1 matrices at this 
specialization, so its determinant has an (n — A:) -fold root at 
q = p k . Likewise T(p,p~ k ) also has rank k and y(qp k )"~ k 
also divides T l . All four of the determinants have this behav- 
ior. In each case, the singular values of q can be read from the 
product formulas for the determinants. The only detail that 
changes is the form of each rank 1 term, which is summarized 
in Table | 

Finally the ^-independent factor R(p) can be found by ex- 
amining the coefficient of the leading power of q, or equiva- 
lently, taking the limit q — » °°. For example 



1 



T i(p,q)ij 

q("+j-i)/2 y(p"+j-') 



i ( " +i+j y 2 C(x,y) 



as q 



• with 



yj = P 



n+j 
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Matrix 



Rank 1 terms 



Extra q value 



p-lipt(n+j) 
p-Uptj 



(p-t> - pf')[ p l{n+j) - p -Z(n+j) ) 



T 4 (p- a -p a ){p- 



In other words, after rescaling rows and columns, TJp,q,r) 
has a block matrix limit: 



lim r'T 5 (p,q,r) 






T x (p,q) 


-T Y {p,q) T 






Table 2: Details of factor exhaustion for Theorem F^, 



In this case R(p) is given by detCj in Theorem |l5|. This hap- 
pens in each case, although for the matrix T 2 it is slightly more 
convenient to specialize to q = 1 . The best extra value of q in 
all four cases is given in Table ^. □ 

Finally we evaluate two three-variable Pfaffians which are 
like the determinants in Theorem [To] 

Theorem 17. For i,j < 2n, let 



T 6 (p,q,r) u = r( P i+i )r( P n 



riq i+i ) 7(^'-0 



Y(pJ+<) y{p]- 1 ) 
y{r j+i ) y(V w ) 
y(pj+i) _ Yip'' 1 ) 



when i ^ j and 



Then 



Pf7 5 = 



P«6 = 



T s(P^^)i.i = ° 
T 6 iP^, r )i,i = 0. 

U, < jYiP M ) 4 U l j7iqP i - i )7(rp J - 1 ) 

n,< ; < 2 „ Y&- 1 ) 
n (<; <2„ r(p J -')n Lj < 2n+l r(qp J -')r(rp J -') 

2\i 

T\i<j<2ny(P i+i ) 



Proof. Factor exhaustion in both q and r. If < k < n, then 

T 5 ip,p\r) id 



E 



<-<e<± 



is, as written, a sum of 2k rank 1 matrices. Therefore the 
Pfaffian, whose square is the determinant, it is divisible by 
liqp )" • Th e same argument applies to T 5 {p,p~ k ,r). It 
also applies to T 5 (p,q,p ±k ) since T 5 is symmetric in q and r. 

This determines Pf T 5 up to a factor R(j>) depending only on 
p. This factor can be determined by taking the limit r — > °°: 



lim 



T s(P>Q> r )i.j 

r {\i-n-\\+\j-n-\\)/2 



T \{P^)j,i-n j< n < 1 





(The bullet • is the exponent above that is different for dif- 
ferent rows and columns.) This establishes that the leading 
coefficient of Pf T 5 (/?, q, r) as a polynomial in r is 

(-l)(5)detli(p,f), 

which in turn determines Rip). 

The Pfaffian Pf T 6 is argued the same way. To find the fac- 
tor Rip) which is independent of q and r, we take the limit 
r,q — > °°. In this limit Pf T 6 reduces to a special case of PfC 4 
in Theorem |l5| □ 

Remark. Several other specializations of the determinants in 
Theorem |l6| and the Pfaffian in Theorem [l7| are special cases 
of Theorem |l5] and other determinants and Pfaffians such as 
these @: 

detjxf 1 } det^jtf- 1 )} 

For example, T { iq 2 ,q) is also a Cauchy double alternant, 
while Tj ip,p") is the product of two (rescaled) Vandermonde 
matrices. Any of these intersections may be used to determine 
the ^-independent factor in the factor exhaustion method. 
There are also other round determinants like the ones in The- 
orem p~q which we do not need, for example 



det(4^ 



Yiq" 



+7+1-1 



Y(pn+j+i-l) 



These examples suggest the following more general problem: 
Let M be an n x n matrix such that M i ■ is a rational polynomial 
in a fixed number of variables, such as p, q, and r, and in 
exponentials of them such as p', q->, and r". When is detM 
round? What if M. • is a rational polynomial in variables such 
as x i andy ; ? 



7. ENUMERATIONS AND DIVISIBILITIES 

In this section we relate the quantities appearing in the other 
sections to prove the results in Section 1 . 
Let 



1 = (1,1,1,1,... ,1) 
x = a 2 + 2 + a 2 
y = aiba)/oiba) 
z = a(ca)/(7(cfl). 



(1) 



Then most of the generating functions in Section 1 can be 
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expressed in terms of the partition functions in Section g, 

Z(n; 1,1) 



A(n;x) = 
A HT (2w;;t,±l) = 
A v (2n;x,y) = 
A m (4n;x,y,z) = 
A QT (4n;x) = 
A Q (2n;x) = 



A vo (8n;x,z) = 



a (a)" 2 -"a (a 2 )' 1 

Zh T (w;T,T) 
a (a) 2 " 2 -"a {a 2 )" 

(7(fl) 2 " 2 -' I a(fl 2 ) n (7(M)" 
z uu(»;lJ) 

a(fl) 4 " 2 -' , (7(a 2 )"(7(/7fl)"(7(cfl)" 

Z QT (n;T) 



(2) 



(7(a) 4 " 2 -"a(a 2 )" 

z («;T) 
a {a) 2 " 2 - 2 " a (a 2 )" 

Z uo (»;T,T) 

C7(fl) 8 " 2 - 3 "(7(fl 2 )"(7(cfl)"/7 2 " ' 



by the definition of the partition functions and the correspon- 
dence between square ice and alternating-sign matrices. The 
generating function A QO (4n;x,y) requires a slightly different 
change of parameters: if y = b 2 /c 2 , then 



A OQ (4n;x,y) 



ZooM) 



o(a) 4 " - 3 "a(a 2 ) n c 2 



The generating function A u (n;x,y) is a polynomial of degree 
n in y and it is easy to show that the leading and trailing coef- 
ficients count VSASMs, so we can say that 

A v (2n + l;x) = A v (n;x,0) =A v (n;x,°°), 

where by abuse of notation, if P(x) is a polynomial (or a ratio- 
nal function), P{°°) denotes the top-degree coefficient. Like- 
wise A uls (n;x,y,z) has bidegree (n,n) in y and z and the corner 
coefficients count VHSASMs and VHPASMs: 

A VHP (4« + 2;x) =A ULJ (n;x,0,°°) =A m (n;x,oo,0) 
A VH (4«+l;x) =A uu (7i;jc,oo,e») 
A VH (4« + 3;x) =A uu (n;x,0,Q). 

We can reverse these relations by defining 

UijaiXifj) (detM± ) 



Zt¥\n;x,y) 



Z^An;x,y) = 



Tlij « fcyj) «(*,-y/) (detM uy ) 

Ui<j a{xpc } )a{y$j) Ui<j o{xp } )a(yp } ) 



7 ( 2)( A , n,<,< 2n o:(x.x ; .)(PfM 00 ) 

n i<; -< 2 „«(^)a(^)(PfM«) 



Z«(«;x) 



n ;<i < 2 „ n,< 7 <2„ ^(^) 



and 

A( 2 )( 2 „;x,±l) 



HT 

A^ v (4n;x,y,z) = o{a) 



:(i( fl r" 2 z±f)(2„;T,T; 



n—2rP~ 



Oiba 



-C7(ca)-"Z( 2 )(4«;T,T) 



A«(4«;x) = (7(«)"- 2 " 2 ZW(4»;T) 
A( ) 2 )(4„;x ; y)=c- 2 ''a(«) 2 "- 2 " 2 z( 2 )(^ 
A[;)(8«;x,z) = (T( fl r- 4 " 2 z( J 1 ( )(8«;T,! 

(2) 

uo 



A. (2 A(8«;x,z 



= C7( fl ) n - 4 " 2 a(cfl)-"z( 2 ) (8n;T,T) 



using the same correspondence between a, b, and c with x, y, 
and z (which slightly different in the case of OOSASMs). The 
observation that all of these quantities must be polynomials 
establishes the factorizations of A HT , A 00 , Aq T , A UU; and A uo 
in Theorem Q 



For vectors x and y, let 

(x,y) = (x v . 



x n,yi i 



,y nj 



denote their concatenation, and let exponentiation of vectors 
denote coordinate-wise exponentiation: 



— {xi,x 2 , ■ ■ ■ ,x k n ). 
Then the matrices 

M{2n;{xX V ),{y,T 1 )) 
M m {2n;{x,x- l ),{y,T 1 )) 

commutes with the permutation matrix 




where /„ is the n x n identity matrix. Similarly the matrices 

M(2« + l;(x,l,^ 1 ),(y,l,y- 1 )) 
M HT (2« + l;(x,l,x- 1 ),(y,l,r 1 )) 

commute with 





(° 





In 


p = 





1 







\I n 









If x has length 2n and b 2 = —c 2 , then 
M 00 (4n;(x,x- 1 )) 

commutes with 

P = 

In each case we can restrict decompose M, M m , and M QQ 
into blocks corresponding to the eigenspaces of P. The block 
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with eigenvalue — 1 is in the three cases equal to M u and pro- 
portional to M ULJ (with b — c — i) and M^, This results in 

the factorizations of A(n;x), A|^|(n;x, 1), and A^(An;x, — 1) 
in Theorem || Another factorization in Theorem |] is that 
of A^,(2n;x 1 — 1). To establish this, observe that the matrix 
M ¥n (2n;x,x) is antisymmetric, and that it is proportional to 
M^,(2n;x). The latter matrix is employed for its Pfaffian, 

while the former for its determinant, which is the square of 
the Pfaffian. 

The final case of Theorem ^ is the relation between 
A y (2n + l;x) and A v (2n;x). This relation is established by 
observing that b appears only in the normalization factor for 
Z u («;x) and not in the matrix M u («;jc); the only step is to 
change variables from b to y. 



Enumeration 


Parameters 




Section p| 


Section p| 


A(n;l) 


a = 






M(q,q(n)) 


tU? 3 ,?) 


A(n;2) 


a = 






M{x,y) 


C,(x 2 ,^) 


A{n;3) 


a = 


a> 6 




M(q,q) 


T 2 (q 3 ,q) 


A HT (2n;l,l) 


a = 


(ft, 




M+ T {q,q{n)) 


7\(</V) 


A HT (2rc;2,l) 


a 






M+ T (q,q) 


T 2 {q\q) 


A v (2n+l;l) 


a — 






M v {q,q(n)) 


T 3 (q 3 ,q) 


A y (2n+1;2) 


a — 






M v (x,y) 


C 2 (x 2 ,f) 


A v (2n + 1;3) 


a — 






M v (q,q) 


T 4 (q\q) 


(4*1,1,1) 


a = 


£U 3 ,i = 


c = ffl 4 


M m (q,q{n)) 


T 3 (q\q 2 ) 


A(9 (4*2, 1,1) 


a = 


fc = c = 


0) 4 




T 4 (q 2 ,q) 


A< 2 > p (4«;l) 


a = 


b = co 3 


c = a 




T 3 (q\q) 


Ag(4«;i) 


a = 






M§(?) 


T 5 (q 3 ,q,q) 


A$(4«;2) 


a = 








T 5 (q\q 2 ,q) 


Ag(4«;3) 


a — 








T 5 (q 6 ,q\q 2 ) 


Ag(4«;l) 


a — 








T 5 (q\q 2 ,q) 


Ag(4«;2) 


a — 


ffl 4 






C 3 (x 2 ) 


A Q (2n;l) 


a — 


ft), 








AW(Sn;l) 


a — 


a> 3 ,c = 


ffl 4 




T 6 (q 3 ,q,q) 


A< 2 >;1,1) 


a = 


tt) 3 ,c = 


ffl 4 




T 6 (q 3 ,q 2 ,q) 



Table 3: Specializations of partition function determinants and Pfaf- 
fians. 



Finally we establish the round enumerations in Theorems [| 
||, and [5| We review the argument from Reference ^ for 
A(n), which is an illustrative case. Let co n = exp(7ri/n), where 
i 2 = — 1. Equation ^implies the following correspondence be- 
tween x and a: 



(On 



a = (0 A 



03, 



x=\ 
x = 2 
x = 3. 



For any of these values of x or a, we would like to evaluate 
the partition function Z(«; 1, 1) to find A(n;x) by equation ^[ 
Unfortunately the matrix M(n; 1, 1) is singular. So instead we 



will find its determinant along a curve of parameters that in- 
cludes (1,1). More precisely, let 



and 



Then 



?W = (^+')/2,^+2)/2 7 ... ; ^+«)/2 ) 



q = q(Q)- 



lim<7(A:) = T, 



and if we assume equation [T], 

M{n;q{0),q{k)) = 



1 



-q k +}-'+X-2-q'-}- k ' 

If we further set x = 1 and A: = n, then 

M(«;^(0)^(«)) = -r i ( ? 3 , ? ) 

has a round determinant by Theorem |l6[ Computing A(n;x) 
is then a routine but tedious simplification of round prod- 
ucts. The argument for most of the other enumerations is 
the same, except that the curve of parameters is (q,q(n)) for 
1 -enumeration in the determinant cases, (q,q) for 2- and 3- 
enumeration in the determinant cases, and q in the Pfaffian 
cases. Each of the matrices is then proportional to some ma- 
trix T i from Theorem 16 or |l7| The determinants and Pfaffian 
for three of the 2-enumerations are round without specializing 
the parameters and instead reduce to a matrix C ( from Theo- 
rem [la. These variations of the argument are summarized in 
TableH. 



8. DISCUSSION 

Even though each section of this article considers many 
types of alternating-sign matrices or determinants in paral- 
lel, the work of enumeratingsymmetry classes of ASMs is 
far from finished. Robbins [ pip conjectures formulas for the 
number of VHSASMs and for the number of odd-order HT- 
SASMs, QTSASMs, and DASASMs in addition to the enu- 
merations that we have proven. Theorem [To] yields a determi- 
nant formula for the number of VHSASMs, obtained from the 
more general partition function Z vv (n;x,y) by setting a = a> 3 
and b — c = flfcf 1 , and all other parameters to 1). In the enu- 



meration of 4n + 1 X An -f 
experiments indicate that 



1 VHSASMs, where b = c 



detMuu^--),^--)) 

is round, but we cannot prove this. In the other case, An + 3 x 
An + 3 VHSASMs, where b = c = a> -1 , we could not even 
find a curve for which the determinant is round. This strange 
behavior of VHSASMs is one illustration that although we 
have put many classes of ASMs under one roof, the house is 
not completely in order. 

For the other three classes conjecturally enumerated by 
Robbins, we could not even find a determinant formula. 
Nonetheless we conjecture: 
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Factor 


;j = l 


n = 2 


n = 3 




l 


x + 2 


x 3 +6x 2 + l3x + 6 


A y (2n;x) 


l 


x + 6 


x 3 + l2x 2 +70x + 60 


A y ^> (4n;x) 


l 


x + 3 


x J + 8x z + 25x + 15 




2 


6x + 4 


20x 3 + 60x 2 + 52x + 8 


AW(4n;x,l,l) 


x + 4 


x 3 +9x 2 +40x+l6 


x 6 + 16x 5 + 125x 4 + 629x 3 + 1036x 2 + 560x + 64 


A®(4n;x) 


l 


x+ l 


x 3 +4x 2 +5x+l 


A VHp( 4 " + 2;X ) 


l 


2x+l 


5x 3 + 12x 2 + 8x+l 


A (2n;x) 


l 


3 


x 2 + 10x+15 


Ag£(4n;x,l) 


2 


2x + 4 


2x 3 +20x 2 +28x + 8 


A«(8n;x) 


l 


x 2 + 5x + 3 


x 6 + 14x 5 + 82x 4 + 210x 3 + 239x 2 + 1 15x + 15 


AW(8n;x) 


2 


2x 2 +20x + 20 


2x 6 + 42x 5 + 420x 4 + 1680x 3 + 2892x 2 + 2040x + 360 


AW(8n;x,l) 


x + 4 


x 4 + l2x 3 +65x 2 + l04x+l6 


x 9 + 24x 8 + 275x 7 + 1966x 6 + 8215x 5 + 19144x 4 



+21777x 3 + 10028x 2 + 1712x + 64 



Table 4: Irreducible x-enumerations 



Question 18. Can DSASMs, DASASMs, TSASMs, and odd- 
order HTSASMs and QTSASMs can be x-enumerated in poly- 
nomial time? 

The polynomials listed in Table Q appear to be generating 
functions of some type, but in most cases there is not even a 
proof that their coefficients are non-negative. (We have more 
data than is shown in the table; the multivariate polynomials 
A^(4n;x,y,z), A^(4n;x,y), and A^(4n;x,y) also appear to 
be non-negative.) Some of them are conjecturally related to 
cyclically symmetric plane partitions pT|]. Indeed they are re- 
lated to each other in strange ways. For example Theorem || 
establishes that if we take x = 1, three of the polynomial se- 
ries (A v , A , and A^ p ) become equal, as if to suggest that 
VSASMs can be .^-enumerated in three different ways! 

Question 19. Do the polynomials in Table ^| x-enumerate 
classes of alternating-sign matrices? 

OSASMs include the set of off-diagonal permutation ma- 
trices, which can be interpreted as the index set for the usual 
combinatorial formula for the Pfaffian. Like ASMs, their 
number is round. These observations, together with the known 
formulas due to Mills, Robbins, and Rumsey p7| ] motivate the 
following question: 

Question 20. Are there formulas for the Pfaffian of a matrix 
involving OSASMs that generalize the determinant formulas 
involving ASMs ? 

Neither any of the enumerations that we establish, nor the 
various equinumerations that they imply, have known bijec- 



tive proofs. Nor is it even known that two equinumerous types 
of ASMs index bases of the same vector space. For example, 
can one find an explicit isomorphism between the vector space 
of formal linear combinations of 2n x 2n OSASMs and the 
vector space of formal linear combinations of 2n + 1 x 2n + 1 
VSASMs? 

Sogo found that Z(«; 1, 1) satisfies the Toda chain (or Toda 
molecule) differential hierarchy Jl7| , [36[ ] . 

Question 21. Ifxandy are set to 1, do the partition functions 
in Theorem [7^ and Table ^ satisfy natural differential hierar- 
chies? 

Many other solutions to the Yang-Baxter equation are 
known [§]. The six- vertex solution corresponds to the Lie 
algebra s 1(2) together with its 2-dimensional representation; 
there are solutions for other simple Lie algebras and their rep- 
resentations. 

Question 22. Do square ice and Izergin-Korepin-type deter- 
minants generalize to other solutions of the Yang-Baxter equa- 
tion ? 

Although our simultaneous treatment of several classes of 
ASMs is not especially short, the argument for any one alone 
is relatively simple. We speculate that the methods of La- 
grange interpolation (used in Section ||) and factor exhaustion 
(the topic of Section |6j) simplify many proofs of product for- 
mulas. I. J. Good's short proof of Dyson's conjecture [ |i"T| | also 
uses Lagrange interpolation. 
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